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Abstract 

The calculation method of the hadron interaction vertex with constant external 
fields in QCD sum rules is discussed. The representation of the polarization operator 
in terms of physical states contributions is considered and the most suitable form 
of this representation is suggested. The estimates of uncertainties in the previous 
calculations of hadron interaction vertices with external fields are given. 



1 Introduction 



The QCD sum rule calculations of hadron interaction vertices with constant external 
fields are widely used for determination of the hadron static properties. In this way the 
proton, neutron and hyperon magnetic moments [1-5] , the axial nucleon and hyperon 
/3-decay couphng constatns [6-8], the nNN, nNA coupling constants [9] and many other 
static hadronic characteristics were found. More recently the same method was used 
to determine the moments of the structure functions [10] and higher twist corrections 
to Gross-Llewellyn and Bjorken sum rules for unpolarized [11] and polarized [12] deep 
inelastic lepton-nucleon scattering. 

I will remind the main features ot the method. The term of quark interaction with 
external field is added to QCD Lagrangian. This field may be a constant electromagnetic 
field in the case of magnetic moments calculations, a fictititious constant axial potential 
in the case of qa and higher twist correction to polarized Bjorken sum rule etc. The 
polarization operator 11 (p^) of currents with the quantum numbers of hadrons in view is 
considered. It is supposed that is negative and —p^ ^ where Rc is the confinement 
radius, and few terms in the operator product expansion (OPE) of 11 (p^) are calculated. 
For quarks the Dirac equation in the external field is written and only the terms linear in 
this field are retained. The term in the polarization operator linear in the external source 
is given by 

n(p2) = J d^xe'P'' < 0\T{r){x),J j{z)d''z,f){0)}\0 > S, (1) 

where S is the external field, rj, fj are the currents with the quantum numbers of the 
hadron, whose interaction vertex with external current j, we would like to determine. 
(E.g. for the case of baryons, rj{x) is proportional to the product of three quark fields, rj ~ 
^abcgagbgc^ where a, b, c are colour indices.) The important ingredient of the calculation is 
the account of induced by the external vield vacuum expectation values in the OPE. In this 
way the linear in the external field term in the polarization operator, proportional to the 
vertex function with zero momentum transfer r(p^) = r(p^,p^,0), 11 = TS is calculated 
by OPE in QCD. On the other side using dispersion relation, n(p^) is represented through 
contributions of hadronic states. Among these the contribution of the lowest state in the 
given channel is of interest and must be separated. If the contributions of exited states 
can be reliably estimated and are small, then by equating of two representations of 11 (p^) 
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the desirable hadron interaction vertex with constant external field can be found. 

In carrying out this program it is essential to represent correctly the hadronic contri- 
butions to n(p^), or what is equivalent, to r(p^), accounting for all possible terms. The 
omission of some terms may result in underestimation of possible errors in the values of 
hadronic coupling constants with external fields and even in completely wrong results. 
Unfortunately, not in all such calculations the general form of hadronic spectra was used 
and the errors due to this, were estimated. The goal of this paper is to discuss the general 
form of hadronic contributions to Il{p^) in the external field (or r(p^)) and to suggest the 
most suitable form for representation of excited states. For few examples the uncertain- 
ties in the determination of hadronic coupling constants with external fields arising from 
contributions of excited states will be estimated. 



2 Dispersion representation of the vertex function 

In order to get the dispersion representation of the polarization operator in the external 
field n(p^) or the vertex function at zero momentum transfer r(p^,p^, 0) it is convenient 
to start from the case, when the momentum transfer q = P2 —pi is not zero, but is small 
and negative. (All what follows refers to the coefficient function at any Lorentz tensor 
structure.) The general double dispersion representation in variables pl,pl of r{pl,pl, q^) 
at fixed q^ <Q has the form (see, e.g. [2]) 

npipi e) = r r , p^X!^'^'^ ds, ds, + p{pi)f{pi e) + p{pi)f{pi q% (2) 

where Pip"^) is a polynomial. (It can be shown that at q^ < there are no anomalous 
thresholds.) For simplicity it is assumed that the currents rj and fj are Hermite conj- 
gate what correspond to the diagonal matrix element over hadronic states. In this case 
p{si,S2,q^) is symmetric in si,S2 and the functions / are the same in the 2^*^ and 3^*^ 
terms in the r.h.s. of (2). The second and the third terms in the r.h.s of (2) plays the 
role of subtraction functions in the double dispersion relation (2). The function /(p^, g^) 
may be represented by one- variable dispersion relation in p^. 

It is clear that the dispersion representation (2) holds also in the limit q 0, P2 ^ 
p\ = p^, where r(p^) = r(p^,p^, 0) is a function of one variable p^. At first sight it seems 
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that one variable dispersion relation can be written for r(p^). Indeed, decomposing the 
denominator in (2), we can write 



r r ,,,^igp> = / / *ii£i) — L_) . (3) 

^0 Jo (S1-P )(S2-P ) J J S1-S2 \S2-P Sl-P / 



In the first (second) term in the r.h.s. of (3) the integration over Si(s2) can be performed 
and the result has the form of one-variable dispersion relation. Such transformation is, 
however, misleading, because, in general, the integrals 



are ultravioletely divergent. This ultraviolet divergence cannot be cured by subtractions 
in one-variable dispersion relation: only the subtractions in the double dispersion repre- 
sentation (2) may be used. It is evident that the procedures, killing the subtraction terms 
and leading to fast converging of dispersion integrals in standard one-variable dispersion 
representations, like Borel transformation in do not help here. 

Let us consider two examples. The first corresponds to the determination of nucleon 
magnetic moments. In this case the current 77 in (1) is the quark current with nucleon 
quantum numbers 77 ~ e'^^'^ffci'ff and the current j in (2) is the electromagnetic current. 
The simple loop diagram - the contribution of unit operator in OPE - is shown in Fig. 
1. It is clear that the spectral function p(si,S2) in (2), corresponding to the diagram of 
Fig. 1 is proportional to 5(si — 52)- The separation of the chirality conserving structure 
results in the statement that the dimension of p is equal to 2 (see [2,3]). So the general 
form of p(si, S2) in (2) is 



where a is a constant. The substitution of (5) into (2) gives for the first term in the r.h.s. 

of (2) at pI = p"^ = p^ 



In this simple example the dispersion representation is reduced to one- variable dispersion 
relation, but with the square of (si — p^) in the denominator. Of course, by integrating 
by parts (6) may be transformed to the standard dispersion representation. However, the 
boundary term, arising at such transformation must be accounted; it does not vanish even 




(4) 



(5) 




(6) 
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after application of Borcl transformation. This means that, even in this simplest case, 
the representation (2) is not equivalent to one-variable dispersion relation. 

The second example corresponds to determination of twist 4 correction to the Bjorken 
sum rule for polarized deep inelastic electron- nucleon scattering [12]. Here the external 
current j in (1) is given by 

= ^qge^^xaG'^.^X'^q, (7) 

where G^j, is the gluonic field strength tensor. An example of the bare loop diagram is 
shown in Fig. 2. In this case, unlike the previous one, p(si, S2) is not proportional to 
S{si — S2). This stems from the fact that in the discontinuity over pf at q ^ and pi ^ pi 
only the left-hand part of the diagram Fig. 2 is touched and the loop integration in the 
right-hand part still persists. For the selected in ref. [12] tensor structure p(si,S2) has 
dimension 4 and is proportional to siS2- We see that in this example the general form of 
dispersion representation (2) must be used in the limit 

q'^0,pl^pl=p'. (8) 

Let us represent r(p^,p^; 0) in terms of contributions of hadronic states, using (2) and 
separating the contribution of the lowest hadronic state in the channels with momentum 
p. Consider the first term in the r.h.s of (2) at = 0,Pi = pI = P^- As is seen from 
Fig. 3 it is convenient to devide the whole integration region in Si, S2 into three domains: 
I) < si < 1^2 < S2 < W^; II) < si < W^, < < 00; < si < 00, 
< S2 < W^; III) < si < 00, < S2 < 00. Adopt the standard in QCD sum 
rule model of hadronic spectrum: the lowest hadronic state plus continuum, starting from 
some threshold W'^. Then in the domain I only the lowest hadronic state h contributes 
and 

p(si,S2) =G'A'5(si -m')5(s2 -m') (9) 

where m is the mass of this state, A is the transition constant of the hadron in the current 
r}. (For /i-baryon < B\fj\0 >= XbVb where vb is the baryon spinor.) G is the coupling 
constant of the hadron with external field, which we would like to determine from the 
sum rule. In the domain III the higher order terms in OPE may be neglected and the 
contribution of hadronic states is with a good accuracy equal to the contribution of the 
bare quark loop (like Figs. 2 or 3) with perturbative corrections. The further application 
of Borel transformation in p^ essentially suppresses this contribution. 
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The consideration of the domain II contribution is the most troublesome and requires 
an additional hypothesis. Assume, using the duality arguments, that in this domain 
also, the contribution of hadronic states is approximately equal to the contribution of the 
bare quark loop. The accuracy of this approximation may be improved by subtraction 
from each strip of domain II of the lowest hadronic state contributions, proportional to 
— m^) or 6{s2 — m'^). Terms of the latter type also persist in the functions f{p\), f{pl) 
in (2). They correspond to the process when the current fj produces the hadron h from 
the vacuum and under the action of the external current j the transition to excited state 
h ^ h* occurs or vice versa (Fig. 4). At pi = pi = P^ these contributions have the form 



with some unknown function b{s). The term (10) will be accounted separately in the r.h.s. 
of (2). I stress that the term (10) must be added to the r.h.s. of (2) independently of 
the form of bare loop contribution p(si, 52)- Even if p(si, S2) = 0, when the OPE for the 
vertex function r(p^,p^, 0) with zero momentum transfer starts from condensate terms - 
the term (10) may persist (the example of such a situation will be given in Sec. 3). (10) 
may be written as 



The integration domain in (12) may be also devided into two parts: < s < W'^ and 
< s < 00. According to our model the contribution of the first part is approximated 
by /i-state contribution, the second one by continuum. These two parts look like the 
contributions of the first and the second terms in the bracket in (11). The first term in 
(11), which after Borel transformation is not suppressed in comparison with the main 
double pole term, arising from (7), as a rule is accounted in the calculations (see, e.g. 
[1-8, 10-12]) as an unknown parameter, determined from the same sum rule. The second 
term in (11), suppressed in comparison with the first one after Borel transformation, is 
traditionally neglected. 

Now we can formulate the recipee, how the sum rule can be written. At the phe- 
nomenological side - the r.h.s of the sum rule - there is the contribution of the lowest 




(10) 




(11) 



The functions /(p^) in (2) can be represented by dispersion relation as 




(12) 
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hadronic state h and the unknown term (11), corresponding to nondiagonal transition 
h ^ h* in the presence of external field: 

+ d^Ks)—-3, -^-3, + (13) 



(p^ — Jw^ s — m? — m? s — j 

(The coefficient a refiects the possibility that in the function / the ratio of terms, pro- 
portional to ijP' — m^)~^ and (s — p^)^^ may differ from 1 as it takes place in (11).) The 
contribution of continuum, corresponding to the bare loop (or also to the higher order 
terms in OPE, if their discountinuity does not vanish at s — > oo) is transferred to the 
l.h.s. of the sum rule. Here it is cancelled by the bare loop contribution from the same 
domain of integration. As a result in the double dispersion representation of the bare loop 
the domain of integration over si, S2 is restricted to < si, S2 < W^. Finally, apply the 
Borel transformation in to both sides of the sum rule. In the l.h.s. - QCD side - the 
contribution of the bare loop has the form 



i dsi / dS2p{si,S2) 
^0 



Si — S2 



= 2Pr d.,/" <i.,«£lliHle-«/"\ (14) 

Jo Jo Si — S2 

where P means the principal value and the symmetry of p{si, S2) was used. The r.h.s. of 
the sum rule is equal; 

A2 



M2 Jw^ s-vn? 



where 



A=r<i,J(£)_ (16) 

Jw^ s — m^ 



Two remarks in connection with eqs. (14), (15) are necessary. If the discontinuity 

p(si,S2) of the bare loop is proportional to 5{si — S2), p(si,S2) = p{si)5{si — S2) like in 
the diagram Fig. 1, then eq. (14) reduces to 

In this case at W'^ ^ the continuum contribution is suppressed exponentially and the 
dependence on the value of continuum threshold W"^ is weak. It, however, p(si,S2) has 
no such form and is a polynomial in si, S2, like in the diagram Fig. 2, then, as can be be 
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seen from (14), the bare loop diagram contribution has a powerhke dependence on W"^. 
In this aspect the QCD sum rule calculation lost a part of its advantages in comparison 
with finite energy sum rules. 

In the case, when the double discontinuity of the bare loop diagram p{si, S2) is propor- 
tional to 5{si — S2), this form will be absent in the radiation correction terms. Here (17) 
is invalid and more general expression (14) must be used. This will result to appearance 
of ln(M^^/ — p^) in the final answer - W'^ plays the role of ultraviolet cut off. 

The necessity to account separately the second term in (15) in the r.h.s. of (15) not 
suppressed after Borel transformation comparing with the first term, was stressed in the 
early treatment of the problem in view [1-3]. In [1-3] it was suggested to separate this 
unknown term from proportional to G term of interest by their different preexponent M ^ 
dependence. This program was successfully realized in the most cases of hadron coupling 
constants with external fields determinations. However, the last term in (15), was omitted 
in all calculations, following the first ones [1-3]. This term is suppressed in comparison 
with the accounted second term in (15) by a factor, less than e~(^^~"*^)/^^. In the most 
cases this factor is of order e^^'^ ~ 1/4. Therefore, if A -C GA^/M^ this term can be safely 
neglected. But at /I ~ GX^/M"^ the nonaccounted term may deteriorate the accuracy of 
the results. 

3 Few examples 

3.1 Quark mass term 

Consider the matrix element 



over the proton state \p >. Here u and d are the fields of u and rf-quarks. This matrix 
element was studied recently [13] by the QCD sum rule technique. On the other side due to 
the Hellman-Feynman theorem [14,15] H is related [13] to the part of the neutron-proton 
mass difference, arising from the quark mass difference /i — ma — m„ 



H —< p\uu — dd\p > /2m, 



'p 



(18) 



(m„ - mpjf, = iiH 



(19) 
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in the linear over /i approximation. Therefore H may be calculated in two ways in the 
QCD sum rule approach: 1) considering three point hadron vertex for interaction with 
external constant uu — dd quark field, as was done in [13]; 2) considering the quark mass 
dependence of proton and neutron mass in the framework of nucleon mass calculation in 
QCD sum rule, as was done in [16,17]. In the latter case the quantity (m„ — rrip)^ - the 
l.h.s. of (19) is calculated. It is evident, that since both calculations are based on the same 
physical ideas and using the same technique, they must be in on-to-one correspondence. 
The comparison of the QCD sides of the sum rules, found in refs. [13] and [17] shows that 
it is indeed the case: the QCD sides of the sum rules for (m„ — mp)^, obtained in [17], 
identically conincide with QCD sides of the sum rules for /iif, found in [13] Q. However, 
the phenomenological sides of the sum rules in refs. [13] and [17] are different. 

The origin of various terms in the phenomenological part of QCD sum rules in the 
calculation, when the quark mass dependence of the nucleon mass was studied, can be 
easily understood. Consider, for example, the chirality conserving sum rule for nucleon 
mass. After Borel transformation the phenomenological part of the sum rule take the 
form 



CO 



A^e-'" + c / e''/'' s^ds, (20) 

where c is some QCD calculated numerical coefficient. We are interested in the dependence 
of (20) on the quark mass difference /z in the linear in /i approximation. So, differentiate 
(20) over /x. We have 

_ A^^e-'"'/^' + dXle-"^"'^'" - c5W\-'^""'"w\ (21) 

where dW"^ is the variation of continuum threshold. Compare (21) with the phe- 
nomenoloigcal side of the sum rule for the vertex iJ, which has the form (15). If the 
last term in (15) is neglected, as was done in [13], then we have 

H_^^mVM^ _ ^^-mVM\ (22) 

The first term in (22) corresponds to the first proportional to 5m?, term in (21) and 
completely agrees with the Feynman-Hellman-theorem. The second term in (22), pro- 
^There was an error in ref. [17] in the coefficient in front of the square of quark condensate term in 
chirality violating structure, instead of 4/3 it must be 2/3. This error was noticed in ref. [f 3]. The results 
of ref. [17] are unaltered. 
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portional to unknown constant A, corresponds to the second term in (21), where 
is also unknown and as well as A is determined from the same sum rule. But nothing 
corresponds to the last term in (21) in eq. (22). This term is proportional to the variation 
of continuum threshold with the quark mass. It is evident that the continuum threshold 
as all hadronic masses is varying with the quark masses. Therefore, the last term in (21) 
generally is non-zero. The contradiction between eqs. (21), (22) disappears if the last 
term in (15) is accounted. The functional dependence of this term is the same as one 
of the last term in (21), since at » JVP the integral in (15) is concentrated at low 
limit 

The requirement of one-to-one correspondence of two sum rules - the sum rule for 
vertex function and the sum rule determining the mass difference, following from the 
Hellman-Feynman-theorem, results in a specific form of the last term in (15), if the stan- 
dard form of hadronic spectrum is used in the sum rules for mass difference. But this is 
not a general case: generally the last term in (15) is arbitrary. 

We see on this example that the last term in (15) must be accounted. This introduces 
an uncertainty in the calculation of three point vertex functions in the QCD sum rule 
approach. 

Consider now another example of the quark mass term contribution - the matrix 
element 



where V is the vector meson, built from qq pair, say ip-meson. The Hellman-Feynman- 
theorem relates to (p — uj {oi tp — p) mass difference 



Due to chiral invariance the OPE for polarization operator in the chirality violating ex- 
ternal field qq starts from the operator mt^qq and the first term in the OPE corresponds 
to the diagram Fig. 5. The dimension of this operator is 4, its contribution vanishes 

at large s and the structure function p(si,S2) in (2) is zero. However, the second term 
in eq. (15) is nonzero. This is evident, if we repeat the same consideration for meson 
mass determination, which were done for nucleon mass and led to (20), (21). Again, the 
requirement that the Hellman-Feynman-theorem holds, leads to the statement that the 



(23) 



(m^ - m^) = TRsH^ 



(24) 
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last term in (15) must be retained. There are no reasons to neglect it, since we are sure 
that the masses of excited states in (/9-channel are heavier than in a;-channel, and, as a 
consequence, the term proportional to SW^ in (21) exists. This example demonstrates the 
necessity to account the term with nondiagonal transitions /i — > /i* in its general form, 
given by the last two terms in eq. (15) even if the bare loop diagram is absent. 

3.2 Nucleon magnetic moments and the axial coupling constants 

9A 

In these cases, [1-8] the bare loop diagram is proportional to 6{si — S2) and the correct 
form of continuum contribution corresponds to (si — p^)^ in the denominator, as in eq. 
(6). This form was used in ref. [2] (see Note Added in Proof) and in refs [6-8]. The 
nondiagonal constant term (17) corresponding to iV — > N* transitions in the external 
field - the second term in (15) - was accounted in the calculations of nucleon and hyperon 
magnetic moments and qa [1-8], but the last term in (15) was neglected. Let us estimate, 
how its account can influence the results. 

In the case of nucleon magnetic moments the constants Ap^^ and Bp „, corresponding 
to the constant A in (15) were found from the sum rules in [2] (see Note Added in Proof). 
Taking these values and assuming that the third term in (15) is suppressed in comparison 
with the second one by the factor cxp[-{W^ - m'^)/M% W'^ = 2.3GeV^ ^ lGeV^ 
we find that the uncertainty in the proton and neutron magnetic moments, arising from 
this source is |A//p| fa |A//„| 0.015, much below the overall uncertainty, estimated in 
[2] - |A/ip,„|/|//p,„| ^ 10%. 

The similar procedure can be performed in the case of isovector nucleon axial coupling 
constant qa determined in [6] and octet axial coupling constant g\ found in [8] . The result 
is the same: the uncertainty, arising due to nonaccounted in [6,8] terms, corresponding to 
the last term in (15) are much less (about 3-4 times) than the extimated overall error. 
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3.3 Twist four correction to Bjorken sum rule for polarized deep 
inelastic lepton-nucleon scattering 

As shown in [18] the twist four correction to Bjorken sum rule for deep inelastic polarized 
lepton-nucleon scattering is expressed through the matrix element over nucleon state of 
the operator 

U, = ugGl^^^^X^u - dgGl^^^-X^d, (25) 

where = {l/2)e^^xaG'\a is the dual gluonic field strength tensor. The calculation of 
the matrix element < p\U^\p >= << U » in the framework of QCD sum rules was 
performed by Balitsky, Brown and Kolesnichenko [12]. The bare loop diagram for this 
case is shown in Fig. 2. In [12] this diagram was calculated, by introducing an ultraviolet 
cut off /i^. It was found that for the chosen Lorentz structure the singular in term 
is proportional to p^ln^(/x^/ — p^)-. Such cut off dependence reflects the fact that the 
spectral function p{si,S2) in eq. (2) is not proportional to 5-function. The logarithmn 
square dependence of r(p^) on the cut-off cannot be removed by Borel transformation. For 
this reasons, in order to obtain physical results the authors of ref. [12] considered various 
values of /i^ in the interval 0| < jj? < IGeV^ and included uncertainties arising from this 
procedure into the error. From the presented above point of view such an approach is not 
legitimate. In this case p(si, S2) is proportional to SiS2'- 

p{si,S2) = bsiS2 (26) 

where 6 is a constant. In the model of hadronic spectrum accepted in Sec. 2, we have after 
Borel transofrmation and using eq. (14): 

Jo Jo Si — S2 



Jo 



s 



(27) 

'0 

Eq. (27) essentially differs from the corresponding expression for bare loop contribution 
in ref. [12]: e.g. the integrand in (27) is positive, while in [12] it is negative in the main 
region of integration. Of course, the QCD sum rule calculation in this case has a serious 
drawback: the dependence of the result on the continuum threshold W'^ is not in the 
form of a small correction of the type cxp(— 14^^/M^) <S 1 at W'^ » M^, but much more 
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strong. This is a direct consequence of high (equal to 5) dimension of the operator ?7^. It 
is clear that the higher the dimension of the considered operator is, the stronger will be 
the dependence on the continuum threshold and less certain the results of QCD sum rules 
calculations. It must be stressed that for operators of high dimensions the loop diagrams 
are in principle nonrenormalizable, the role of excited states in the physical spectrum 
increases and the determination of the lowest state contribution becomes impossible. 

Eq. (27) must be taken instead of the contribution of the bare loop diagram, used 
in [12]. A similar procedure must be also applied in the case of other terms in the sum 
rules [12], containing ultraviolet cut off. The results of such calculation show, however, 
that occasionally the numerical value of <^ f/ ^ obtained in this way does not differ 
essentially form the value found in [12] - it is in the limit of quoted errors. (The same 
refers to the other calculated in [12] matrix element <^V related to the integral from 
the structure function g2{x).). Unfortunately, this fact does not mean that the value of 
<^ U ^ determined in [12] is rehable. The value <S C/ » obtained from the sum rule 
almost entirely comes from the last accounted in [12] term of OPE of dimension 8. (The 
contribution of the bare loop corresponds to the background term - the second term in 
eq. (15).) For this reason unless the next term of OPE will be calculated in this problem, 
there is no confidence in the result. 

4 Conclusions 

It was demonstrated, that in the QCD sum rules determination of hadronic coupling 
constants with external fields the phenomenological sides of the sum rules were not treated 
properly in some cases. The most suitable form of the representation of the physical 
states contributions - eq. (15) is suggested. At the same time in the loop diagrams in 
the QCD sides of the sum rules the integration in the double dispersion representation 
must be performed in the intervals < s S2 < W^, where VT^ is the continuum threshold. 
It is shown that for high (negative in mass) dimension of external field the results of 
coupling constants determination in essential way depend on the value of the continuum 
threshold W'^ and are less certain than other QCD sum rule results. It was stressed 
the necessity the account in the nondiagonal transition terms, corresponding to diagram 
of Fig. 4, not only nonsuppressed by Borel transformation terms, but also the terms. 
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exponentially suppressed by Borel transformation - the last term in eq. (15). The role of 
these nonaccounted in previous calculation terms was estimated for few examples. 
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Figure captions 



Fig. 1. The bare loop diagram, corresonding to the determination of nucleon magnetic 
moments or to the quark mass correction to the baryon masses. The solid lines correspond 
to quark propagators, crosses mean the action of currents 77, 77, the bubble corresponds to 
quark interaction with external field. 

Fig. 2. The bare loop diagram for twist 4 correction to the Bjorken sum rule for deep 
inelastic electron-nucleon scattering. The dashed vertical line corresponds to discontinuity 
over pI at 7^ 

Fig. 3. The integration domains in si,S2 plane. 

Fig. 4. The schematical representation of /i — > h*{h* — > h) transitions in the external 
field. 

Fig. 5. The diagram representing the contribution of the operator TJiqqq to the sum rule 
for quark mass correction to vector meson mass. 
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